Holographic Superfluids and the Dynamics of Symmetry Breaking 
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We explore the far from equilibrium response of a holographic superfiuid using the AdS/CFT 
correspondence. We establish the dynamical phase diagram corresponding to quantum quenches of 
the order parameter source field. We find three distinct regimes of behaviour that are related to the 
spectrum of black hole quasi-normal modes. These correspond to damped oscillations of the order 
parameter, and over-damped approaches to the superfiuid and normal states. The presence of three 
regimes, which includes an emergent dynamical temperature scale, is argued to be a generic feature 
of time-reversal invariant systems that display continuous symmetry breaking. 
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In the last few years there has been a wealth of experi- 
mental activity exploring the non-equilibrium properties 
of quantum many body systems. Recent advances in- 
clude observations of long-lived oscillations in colliding 
one-dimensional (ID) Bose gases [Tj, and the dynamics 
of cold atoms following a quantum quench [2] [3] . Non- 
equilibrium measurements have also been exploited to 
reveal the amplitude mode in 2D superfluids [3], and to 
explore pairing in high temperature superconductors [5] . 
In parallel there has also been significant theoretical work 
on strongly correlated systems, particularly in ID, where 
analytical and numerical progress is possible. This in- 
cludes predictions for the time dependence following a 
quench to a quantum critical point [5] , and the develop- 
ment of numerical techniques [7HH]- Significant progress 
has also been made on the quench dynamics of integrable 
models, including the Ising model, the ID Bose gas, and 
the BCS Hamiltonian; for a recent review see Ref. [10] . 

A notable feature of the integrable dynamics of the 
BCS Hamiltonian, following an abrupt change in the 
pairing interactions, is a regime of undamped persis- 
tent oscillations of the order parameter [TTMT5] ; see also 
Refs. [TH [T5] . This is accompanied by a transition to a 
regime of damped collective oscillations as the strength of 
the quench is increased [16] . These results have been ob- 
tained using time-dependent generalisations of the BCS 
wave function and through exact spectral considerations. 
They are expected to hold in the collisionless regime, 
for timescales shorter than the energy relaxation time 
[T7l [18] . In spite of these achievements, it is challenging 
to see how such results are modified at late times in the 
collision dominated regime. In particular, do the oscilla- 
tions and the transition withstand quantum and thermal 
fluctuations and departures from integrability? Indeed, 
related considerations apply to other integrable systems, 
and generalizing non-equilibrium results to more generic 
situations, including higher dimensions, is a major open 
challenge. 

In this respect, the AdS/CFT correspondence [TTjH2T] 



can offer valuable insights. It recasts certain strongly in- 
teracting quantum systems, which are large N field the- 
ories, in terms of weakly coupled gravitational models in 
at least one dimension higher. This provides access to the 
quantum response functions from the classical solutions 
of the corresponding gravitational equations. Moreover, 
finite temperatures can be incorporated by the inclusion 
of black hole solutions. This real time gravitational ap- 
proach therefore avoids the analytical continuation re- 
quired by the imaginary time path integral. It is thus 
ideally suited to non-equilibrium problems at finite tem- 
perature and in arbitrary dimension |22H27j . The meth- 
ods are very powerful when combined with numerical so- 
lution of the equations of motion, as they allow access 
to the far from equilibrium response over the entire time 
evolution [2"5P5]. 

In this manuscript we will focus on the far from equi- 
librium response of a holographic superfiuid 29, 34, 35 
under a spatially homogeneous and isotropic quench. 
Our primary aim is to reveal the presence of three 
regimes of non-equilibrium response, including a dynam- 
ical transition from under-damped to over-damped col- 
lective oscillations. This transition occurs within the su- 
perfiuid phase, and we argue that this transition will fea- 
ture in other (holographic and non-holographic) time- 
reversal invariant systems that display continuous sym- 
metry breaking. 

Model. — We consider the original action for a holo- 
graphic superfiuid, describing a complex scalar field tp, 
with charge q and mass m, minimally coupled to electro- 
magnetism and gravity in 3+1 dimensions |34l I35j : 
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where F a b = d a Ab — dbA a , D a = d a — iqA a , and the radius 
t parameterizes the inverse curvature of AdS space-time. 
The model is dual to a strongly coupled large N CFT in 
2+1 dimensions, residing on the AdS boundary; for re- 
cent reviews see Refs. [30JI3Z]- This CFT is time-reversal 
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invariant and has a global U(l) symmetry whose con- 
served current, J M , is dual to A a . Likewise, the field ip 
is dual to an operator O in the CFT with scaling dimen- 
sion A fixed by m. For simplicity, we choose q = 2 and 
m 2 = —2/£ 2 with A = 2. As we will discuss more fully 
below, the operator O corresponds to a superfluid order 
parameter for the boundary theory |34l 135) . In the re- 
mainder we will denote 1/(2k 2 ) = C/£ 2 and choose units 
with I = 1. Here, C is a measure of the number of local 
degrees of freedom in the CFT, going as C ~ N 3 ^ 2 at 
large N. As we will see below, the model exhibits global 
U(l) symmetry breaking, which is possible in the 2+1 
D CFT due to large N [35 . We emphasise that other 
holographic superfluid models, including 3+1 D CFTs, 
will exhibit analogous phenomena. 

In order to investigate the dynamics of the model we 
need to specify the space-time coordinates and the met- 
ric. Assuming spatial homogeneity and isotropy of the 
CFT under time evolution, the most general metric is 

ds 2 = z- 2 [-Fdt 2 -2dtdz + S 2 {dx\ + dx 2 2 )} , (2) 

where (t, Xi,X2) are common to both the boundary and 
bulk theories, and the coordinate z specifies the addi- 
tional direction in the bulk space-time. This corresponds 
to ingoing Eddington-Finkelstein coordinates, where the 
asymptotic AdS boundary is located at z — 0; see Fig. [TJ 
Here F(t, z) and S(t, z) depend only on t and z. Likewise 
ip = tp(t, z) and At = A t (t, z), where the spatial compo- 
nents of A a are set to zero. We may thus examine the 
dynamics as a function of t and z alone; see Fig. [I] 

First, let us recall that in equilibrium, the CFT at 
finite temperature and charge density is described by an 
electrically charged static black hole whose properties are 
independent of t. The high temperature unbroken phase 
of the equilibrium CFT is described by an AdS-Reissner- 
Nordstrom (AdS-RN) black hole [3H ES] with 

F = 1 — 2Mz 3 + (p 2 /4)z 4 , S=l, A t = fj, - pz, (3) 

and "0 = 0. Here p and p are the chemical potential and 
the charge density of the dual CFT; the latter is related 
to the expectation value of the current density operator 
via the relation (J t ) = p/2n 2 . Likewise, the mass M is 
proportional to the energy density of the CFT. Similarly, 
the Hawking temperature, Ta(p,M), corresponds to the 
temperature T of the CFT. At T c w 0.090// the AdS-RN 
black hole becomes unstable and the CFT is described 
by a new family of black hole solutions with ip ^ [35] . 
Asymptotically close to the AdS boundary ip has the co- 
ordinate expansion ip = zipi + z 2 ip>2 + ■•• with ipi = 0. 
Analogous to the identifications following Eq. ([3]), the 
AdS/CFT correspondence allows one to identify ip\ as 
a source for the superfluid order parameter in the dual 
CFT, and ip2 as the expectation value, (O) = tp2 /2k 2 . 
Hence, these black holes with ip ^ describe a super- 
fluid phase in which the global U(l) symmetry is spon- 
taneously broken. 



Gaussian Quantum Quench. — We now analyse the far 
from equilibrium dynamics of the dual CFT, at finite 
temperature and charge density, by numerically con- 
structing time-dependent black hole solutions, where 
F, S,At,ip depend on both z and t. Near the AdS bound- 
ary at z = the latter have the asymptotic expansion 

ip = zip 1 {t)+z 2 ip 2 {t) + ..., A t = p(t)-zp(t) + ... (4) 

The holographic renormalization group allows one to es- 
tablish the time-dependent correspondence 

Wtw) - — ^2 — ' \ U W/ ^2 ' 

(5) 

in our space-time coordinates and gauge. We take as 
our initial state a superfluid corresponding to a black 
hole with ip |3S]. We set the initial temperature 
to Tj = 0.5T C for numerical convenience; as we will see, 
similar results are also expected for other values of TJ. We 
then apply a quench of the source field tpi(t), conjugate 
to (0(t)). Specifically, we apply a Gaussian- type quench, 
centred on t = 0, by imposing 

^(t) = 5e-^\ (6) 

where S and f characterise the quench strength and time- 
scale respectively. The chemical potential of the initial 
state, pi, sets the scale for the resulting dynamics and 
we use it to define dimensionless S = S/pi and r = pi?. 
For definiteness, we set t = 0.5 and will vary S. We 
track the dynamics by solving the equations of motion 
of ([T]) numerically. As discussed in the Supplementary 
Material we choose a gauge for p(t) which keeps the ini- 
tial and final charge densities the same while the quench 
injects energy into the system. Then p{t) interpolates 
from the initial value, pi, to a final chemical potential, 
Pf. We find similar results for other values of r, and 
also for quenches which do not preserve the equality of 
initial and final charge density. As we shall see later, our 
quench is abrupt compared with the emergent relaxation 
time-scale. Further details on solving the equations of 
motion, which are non-linear PDEs in t, z, are provided 
in the Supplementary Material. Here we focus on the 
main results, summarised in Figs[T]|3] 

Dynamical Phase Diagram. — In Fig. [2] we show the 
dynamical phase diagram as a function of the quench 
strength, 5. It displays three regimes of late-time be- 
haviour whose asymptotics are governed by the gauge- 
invariant equation 

KO(t))|*KO>/ + .Ae- fa ' t | 1 (7) 

where (O) / is the final order parameter, A is an ampli- 
tude pre-factor, and a; is a complex frequency in the lower 
half-plane. In region III it displays exponential decay to- 
wards a vanishing final order parameter = 0, so 
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FIG. 1. Schematic representation of the space-time coordi- 
nates Q. We show data for the time evolution of the real 
part of the scalar field i[>(t, z) following a Gaussian quench at 
t — with 5 — 0.15, from a superfluid black hole initial state 
as t — > — oo with Ti/T c — 0.5. The behaviour near the AdS 
boundary at z — is used to extract the dynamics of the 
superfluid order parameter (0(t)) shown in Figs 2[3 



that for large quench strengths we exit the initial super- 
fluid phase completely. In contrast, in region II it exhibits 
non-oscillatory exponential decay with Re(w) = to- 
wards |(0) /| 7^ 0. As we shall see later, this corresponds 
to the presence of a gapped "amplitude" mode and a 
gapless "phase" mode in the superfluid phase. However, 
in region I it exhibits exponentially damped oscillations 
with Re(cj) ^ towards 7^ 0, so that for smaller 

quench strengths there is another regime of dynamics. 
For the parameters used in Fig. [2j the transition from 
I to II occurs at a critical quench strength 5* ~ 0.14, 
whilst the transition from II to III occurs at S c sa 0.21. 

The behaviour shown in Fig. [2] is reminiscent of the 
dynamical phase diagram for a BCS superconductor 
[16], despite the fact that the holographic superfluid is 
strongly coupled, and that the effects of thermal damp- 
ing are incorporated. Indeed, the persistent oscillations 
of the integrable BCS Hamiltonian are replaced here by 
an under-damped approach towards 0, whilst 

the power-law damped BCS oscillations are replaced by 
an exponentially damped approach. The transition at 
5* provides a finite temperature and collision dominated 
analogue of the collisionless Landau damping transition 

Emergent Temperature Scale. — We can gain further 
insight by considering the phase diagram as a function 
of the final temperature, Tf, corresponding to the equi- 
librium temperature of the final state black hole. In 
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FIG. 2. (a) Dynamical phase diagram of the holographic su- 
perfluid showing the final order parameter, at late 
times. We start in the superfluid state with Ti — 0.5T C , and 
monitor the time evolution with increasing quench strength 
S. The dynamics exhibits three regimes of behaviour; for the 
chosen parameters the transitions between these regimes oc- 
cur at 5* ~ 0.14 and S c « 0.21. (b) In region I we observe 
damped oscillations towards 7^ 0. (c) In region II we 
find a non-oscillatory approach towards 7^ 0. (d) In 
region III we find a non-oscillatory decay towards \ {0) j\ — 0. 
The dashed lines in (b), (c), and (d) correspond to the dom- 
inant quasi-normal modes of the final state black holes for 
temperatures Tf /T c = 0.73,0.95, 1.48 respectively. 



Fig. we plot Tf versus <5, showing that stronger 

quenches lead to greater final temperatures, consistent 
with the notion that the quench leads to heating. Using 
this relationship we may re-plot the data in Fig. [2] as a 
function of Tf] see Fig. j3^b) . The data collapse on to the 
equilibrium phase diagram of the holographic superfluid 
[35], as indicated by the solid line. The transition from 
II to III is associated with increasing Tf above T c . How- 
ever, Fig. [3jb) contains more information than the equi- 
librium phase diagram; there is an emergent dynamical 
temperature scale , associated with (5* , where the time 
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FIG. 3. (a) Final state temperature T/ versus the quench 
strength S. We use the same initial parameters as in Fig. [2] 
(b) Using the relationship in (a) the dynamical boundary in 
Fig. |2| collapses on to the equilibrium phase diagram of the 
superfluid 34, 35 as shown by the solid line. The dynamical 
transition at S t ~ 0.14 in Fig. [2] occurs within the superfluid 
phase at a temperature T, ~ 0.81T C . 



evolution changes from under-damped to over-damped. 

Quasi-Normal Modes. — In order to gain insight into 
the three regimes of collective dynamics, the temperature 
scale T* , and to expose its broader relevance, we turn our 
attention to the late time asymptotics in Figs 2]|3 As 
t — ^ oo, the dynamics is described by the quasi-normal 



modes (QNMs) of the late time black holes. Each QNM 
describes an approach to equilibrium in linear pertur- 
bation theory with time dependence e~ tut . Those that 
dominate the late time dynamics have complex frequency 
uj with the greatest imaginary part and give rise to the 
behaviour in Eq. Q. As outlined in the Supplementary 
Material, we have calculated the homogeneous isotropic 
QNMs both for the AdS-RN black holes (jS), and for the 
superfluid black holes of Ref. [35] . The trajectories of the 
dominant QNMs in the complex uj plane are depicted in 
Fig. [4] 

Typically, the real parts of the dominant QNM fre- 
quencies correspond to oscillations, and the imaginary 
parts to damping. However, as shown in Fig. |4]jc), for 
temperatures above T c , the QNMs for the AdS-RN black 
hole have two complex frequencies that are closest to the 
real axis. Nonetheless, substitution into Eq. ^ with 
(Of) = yields the damped non-oscillatory behaviour 
found in region III of Fig. [2j As the temperature is low- 
ered, these dominant poles migrate upwards in the com- 
plex uj plane and at the superfluid transition tempera- 
ture, T c , they coincide at co — 0. This corresponds to the 
onset of spontaneous U(l) symmetry breaking with the 
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FIG. 4. Evolution of the QNM frequencies with temperature, 
(a) T = 0.71T C . (b) T = 0.95T C . (c) T = 1.68T C . Time 
reversal invariance corresponds to cj — > —uj*. (d) and (e) 
show the imaginary and real parts of the dominant QNMs. 
i.e. the QNMs closest to the real axis. The results show three 
regimes of dynamics, in quantitative agreement with Fig. [2] 



appearance of a Goldstone mode. Below T c , one of these 
modes, the "amplitude" mode, travels down the imagi- 
nary axis, consistent with time-reversal invariance under 
uj — > — oj*, whilst the Goldstone "phase" mode remains 
pinned at uj = 0; see Fig. [4jb) . The amplitude mode de- 
scribes the damped approach to a finite order parameter 
as shown in region II of Fig. [2j the Goldstone mode does 
not affect the dynamics in the homogeneous and isotropic 
context, although it does lead to a hydrodynamic mode at 
non-zero spatial momentum. As the temperature is low- 
ered, the subdominant poles also ascend in the complex 
plane. The dynamical transition temperature corre- 
sponds to the temperature at which the damping rate 
of the descending amplitude mode coincides with that of 
the ascending subdominant complex poles. For the cho- 
sen parameters this occurs at w 0.81T C , in agreement 
with the nonlinear analysis. For temperatures below T*, 
the previously subdominant poles now become dominant, 
as shown in Fig. |4|a). The dynamics corresponds to a 
damped oscillatory approach to a finite order parameter 
as found in region I. In addition to this change in dynam- 
ics at , one may also extract the variation of the emer- 
gent timescales as a function of temperature. As shown 
in Figs.Qd) and (e), there are three regimes. Moreover, 
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the extracted timescales are in quantitative agreement 
with the late time behaviour of the nonlinear analysis, 
as indicated by the dashed lines in Figs. [2jb)-(d). It is 
readily seen that the linear response analysis provides an 
excellent description over a broad time interval. 



Dynamics of Symmetry Breaking. — We now argue 
that the main results on the late time behaviour of the 
quenched holographic superfluids, captured in Figs. [2]|4j 
have a surprisingly universal applicability. To see this 
we first recall, using the holographic dictionary, that the 
location of the QNMs of the black holes, presented in 
Fig. |4j correspond to the location of poles of the retarded 
Green's function for the operator O in the dual theory 
[55] , That is, the results can be reformulated in terms 
of the dual CFT itself, without reference to black hole 
solutions. Thus, the late time behaviour of the system 
is governed by the poles of the retarded Green's function 
that are closest to the real axis, exactly as we have ar- 
gued above. In particular, the appearance of the scale 
T» and the three types of behaviour, labelled I, II and 
III in Figs. [2] and [3J are a consequence of the pole struc- 
ture presented in Fig.J3] Now the key point is that the 
pole structure in Fig. |4jis the generic behaviour one ex- 
pects for an arbitrary homogeneous and isotropic system 
with time-reversal invariance under uj — > — uj* , which can 
spontaneously break a continuous global symmetry in- 
cluding the presence of the Goldstone mode at the ori- 
gin. The precise value of T* will, of course, depend on 
the details of the system, but it will always be given by 
the temperature at which the value of Im(w) for the pole 
on the imaginary axis, and those poles off the imaginary 
axis and closest to the real axis, coincide. In general we 
will have < X"* < T c and hence all three regions I, II 
and III; for some systems it may be that T* = and 
hence region I would be absent. For a local symmetry we 
expect the same phenomenology; the Goldstone mode is 
replaced by the longitudinal mode of the massive vector 
which, together with the amplitude mode, will be con- 
strained to be on the lower imaginary axis. It may be 
possible to compute the pole structure in other models; 
see for example [3j5] for a recent calculation of the spec- 
tral properties of the O(N) model at zero temperature. 
Finally, it would be very interesting to experimentally 
observe the behaviour we are predicting in, for example, 
a system of cold atoms. 
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SUPPLEMENTARY MATERIAL 

Nonlinear Dynamics. — To study the response to the 
quench given in Eq. ^ we solve the equations of mo- 
tion for Eq. numerically, using the metric described 
near ([2| (see also [25]). The asymptotic boundary is lo- 
cated at z — and writing ip(t,z) — z[ipi(t) + ip(t, z)], 
a(t, z) = fi{t)+a(t, z), F{t, z) = l + z 2 [-\ifj 1 \ 2 /2 + F(t, z)] 
and S(t, z) = 1 + z 2 [—\ipi\ 2 /4 + S(t, z)] we can choose a 
gauge where ip ~ ip2{i)z, a ~ p(t)z, and F, S also van- 
ish linearly as z — > 0. Notice that the residual coordi- 
nate freedom l/z — > 1/z + f(t) is fixed by these asymp- 
totics. In contrast to [25], we adopt a gauge for fi so 
that Im(V>2 — Dipi) = 0, where D = dt — 2ifx. One can 
show that the equations imply the boundary charge con- 
servation equation, given by p = — 4Im['i/'*('02 — Dipi)], 
and also a (sourced) energy conservation equation. Our 
choice of quench, with tpi real as in Eq. ©, ensures that 
the initial and final charge densities are the same. Specif- 
ically, this can be seen by considering Eq (5) and observ- 
ing that we have p = and also fi = at t = ±oo. 
With S, F, A t and the complex scalar ip we have 5 real 
quantities to evolve as functions of the coordinates t, z. 
There are 8 (real) equations of motion from the met- 
ric, vector and scalar equations, 5 of which we use as 
evolution equations for the variables ip,a,F,S, and the 
remaining 3 are constraints. The 5 equations have prin- 
cipal parts, d 2 z ip — \Fd 2 ^) = 0, d 2 z a = 0, df z S = 0, 
and Fd 2 S + Sd^F = 0, where the latter is elliptic. Pro- 
vided the initial data satisfies the 3 constraints then two 
of them, with principal parts d 2 S = and d 2 a = 0, are 
automatically satisfied at later times. The one remain- 
ing constraint, with principal part d 2 S = 0, must be 
imposed at the boundary in addition to the 5 equations 
of motion in the bulk. Physically, this corresponds to 
(sourced) boundary stress energy conservation, and pro- 
vides the additional data required for the elliptic evolu- 
tion equation. Note that (sourced) current conservation 
results from the evolution equations and is not imposed 
separately. 



To solve the equations of motion numerically we use 
a Chebyshev pseudo spectral representation in z £ [0, 1], 
where the initial horizon is located at z = 1, and use 
an implicit Crank-Nicholson finite difference scheme in 
t. Given a slice at constant t we advance to the next 
slice by solving for ip, a, F, S from the 5 evolution equa- 
tions together with the boundary constraint and gauge 
condition for /x; note that we are not required to impose 
any boundary condition at the innermost point z = 1 
since this is inside the event horizon. The method out- 
lined above is very robust, and by virtue of the spectral 
representation in z it allows very accurate extraction of 
boundary quantities. We find that relatively modest grid 
sizes with 20 points in z already give reliable results. The 
data presented in this paper is for 40 points, where con- 
vergence testing indicates that the errors will be less than 
percent level in all the plotted quantities. We construct 
static superfluid black hole solutions of [35] for the initial 
data by solving the usual ordinary differential equation 
shooting problem in our gauge. We then find that we 
can stably evolve up to times t ~ 60 p~ x , to extract the 
data presented here. If one runs for too long then we 
sometimes encounter the black hole singularity. In order 
to evolve further one must implement singularity exci- 
sion which we leave for future work. Our present results 
are in full quantitative agreement with the QNM analysis 
described in the text and below. 

Quasi-Normal Modes (QNMs). — The QNMs of the 
equilibrium black holes of Eq. ([T]) are linearised pertur- 
bations with ingoing boundary conditions at the black 
hole horizon and normalisable boundary conditions at 
the asymptotic AdS boundary; for a review see Ref . [3S] . 
For the purposes of this paper we only consider the zero 
momentum sector of the QNM spectra. Furthermore, to 
determine the late dynamics of the charged scalar field, 
ip, we only analyse the sector involving ip. For the AdS- 
RN black hole ^ we find a second order ODE for ip. For 
the superfluid black hole of [3S] which has ip ^ 0, we can 
use the gauge freedom, arising from diffcomorphisms and 
local U(l) transformations of the background, to reduce 
the problem to two second order ODEs for two gauge in- 
variant variables involving ip, A and the metric. In both 
cases we then use Chebyshev pseudo-spectral differencing 
to cast the linear perturbation equations into the form, 
M(lj; A)v = 0. The matrix M depends on the complex 
frequency uj and the background parameters A, such as 
the temperature and chemical potential. The vector v 
consists of the two gauge-invariant variables evaluated at 
the grid points. The QNM frequencies are then deter- 
mined by the condition |M(w; A)| =0. 



